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Abstract 

A  fuel  cell’s  output  power  depends  nonlinearly  on  the  applied  current  or  voltage,  and  there  exists  a  unique  maximum  power  point  (MPP). 
This  paper  reports  a  first  attempt  to  trace  MPPs  by  an  extremum  seeking  controller.  The  locus  of  MPPs  varies  nonlinearly  with  the  unpredictable 
variations  in  the  fuel  cell’s  operation  conditions.  Thus,  a  maximum  power  point  tracking  (MPPT)  controller  is  needed  to  continuously  deliver  the 
highest  possible  power  to  the  load  when  variations  in  operation  conditions  occur.  A  two-loop  cascade  controller  with  an  intermediate  converter  is 
designed  to  operate  fuel  cell  power  plants  at  their  MPPs.  The  outer  loop  uses  an  adaptive  extremum  seeking  algorithm  to  estimate  the  real-time 
MPP,  and  then  gives  the  estimated  value  to  the  inner  loop  as  the  set-point,  at  which  the  inner  loop  forces  the  fuel  cell  to  operate.  The  proposed 
MPPT  control  system  provides  a  simple  and  robust  control  law  that  can  keep  the  fuel  cell  working  at  MPPs  in  real  time.  Simulation  shows  that  this 
control  approach  can  yield  satisfactory  results  in  terms  of  robustness  toward  variations  in  fuel  cell  operation  conditions. 

©  2007  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

As  a  clean  energy  conversion  technology,  fuel  cells  are  a 
promising  alternative  to  a  wide  variety  of  power  generation 
appliances.  Due  to  its  low  operating  temperature,  high  power 
density  and  fast  startup,  the  proton  exchange  membrane  (PEM) 
fuel  cell  power  plant  is  a  promising  candidate  for  residential  and 
vehicular  applications.  But  fuel  cells  call  for  large  amounts  of 
investment.  To  minimize  the  overall  system  cost,  the  ability  to 
extract  the  maximal  power  from  a  fuel  cell  is  a  crucial  issue  that 
must  be  considered  for  the  optimal  design  of  a  fuel  cell  powered 
system  [1]. 

When  a  fuel  cell  is  directly  connected  to  an  external  load, 
its  output  power  depends  on  both  the  internal  electro-chemical 
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reaction  and  the  external  load  impedance.  The  system’s  oper¬ 
ating  point  is  at  the  intersection  of  the  fuel  cell’s  I-P  curve 
and  the  load  line.  There  exists  a  unique  operating  point,  called 
the  maximum  power  point  (MPP),  at  which  the  fuel  cell  pro¬ 
duces  its  maximum  power,  as  illustrated  in  Fig.  1.  According 
to  the  power  transfer  theory,  the  power  delivered  to  the  load 
is  maximized  when  the  fuel  cell  internal  impedance  equals  the 
load  impedance.  As  either  of  the  two  impedances  may  vary,  the 
chances  are  quite  slim  that  the  MPP  happens  to  be  the  operating 
point.  In  most  cases,  it  is  actually  undesirable  to  operate  at  the 
maximum  power,  because  the  corresponding  fuel  efficiency  is 
at  best  50%.  But  there  are  applications  where  power  density  is 
at  a  premium  as  compared  with  fuel  efficiency,  and  the  maxi¬ 
mum  power  point  tracking  (MPPT)  would  be  beneficial  in  such 
situations. 

Extensive  discussions  on  the  fuel  cell  I-V  characteristic, 
which  can  be  found  in  the  literature  [2-4],  show  that  the  curve 
depends  strongly  on  operating  parameters  in  a  nonlinear  way. 
Mishra  et  al.  did  some  valuable  work  on  the  manipulation  of  the 
fuel  cell  to  attain  the  MPP.  Firstly,  they  picked  out  a  physics- 
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Fig.  1.  Typical  fuel  cell  polarization  and  power  curves. 


based  model  for  fuel  cell  design  and  optimization  [1].  Then, 
they  made  use  of  the  model-based  optimization  methodology  to 
determine  the  MPP  of  the  proton  exchange  membrane  (PEM) 
fuel  cell,  which  is  influenced  by  such  operating  parameters  as 
the  fuel  cell  temperature,  the  anode  pressure,  the  cathode  pres¬ 
sure,  the  relative  anode  humidity,  the  relative  cathode  humidity, 
the  anode  stoichiometry,  the  cathode  stoichiometry,  the  anode 
dry  gas  mole  fraction,  and  the  cathode  dry  gas  mole  fraction  [5]. 
In  essence,  they  took  an  open  loop  approach.  Hence  the  model 
mismatch  caused  by  parametric  uncertainties  and  disturbances 
is  unavoidable.  Also  it  is  very  hard  to  online  determine  whether 
the  fuel  cell  is  working  at  its  MPP  by  measuring  so  many  param¬ 
eters.  Most  importantly,  they  neglected  the  other  side  of  the  coin: 
the  external  load. 

Benziger  et  al.  [6]  noted  that  previous  researchers  had  over¬ 
looked  a  decisive  factor  in  the  determination  of  the  MPP,  i.e.,  the 
external  load.  By  adding  a  variable  load  to  the  base  external  load 
to  achieve  increased  power  output,  they  found  that  changing  the 
external  load  impedance  was  a  more  effective  method  to  alter  the 
power  delivery.  It  is  well-known  that  maximum  power  occurs 
when  the  sum  of  the  external  impedance,  i.e.,  the  variable  load 
impedance  plus  the  base  load  impedance,  is  equal  to  the  fuel  cell 
impedance.  Yet  as  mentioned  above,  the  location  of  the  MPP  in 
the  I-P  plane  is  hard  to  predict.  How  to  control  the  variable  load 
to  track  the  MPP  is  thus  a  rather  challenging  task. 

Up  to  now,  few  researches  on  fuel  cell  MPPT  control 
have  been  reported,  but  a  number  of  methods  have  been  used 
in  photovoltaic  power  applications  [7],  such  as  the  perturba¬ 
tion  and  observe  (P&O)  method,  the  conductance  incremental 
method,  the  parasitic  capacitance  method,  the  only  current 
method,  model-based  methods,  and  artificial  intelligence  meth¬ 
ods.  Among  them  the  P&O  algorithm  is  by  far  the  most 
commonly  used  in  practice  because  of  its  ease  of  implementation 

[8] .  But  due  to  the  limitation  that  P&O  exhibits  erratic  instable 
behavior  under  rapidly  changing  environments,  the  algorithm  is 
unsuitable  for  the  job  of  tracking  the  frequently  moving  MPP 

[9] . 

In  this  paper,  a  maximum  power  point  tracking  (MPPT)  con¬ 
troller  is  designed,  which  can  be  used  to  adjust  the  load  to  the 


source  impedance  so  that  the  equilibrium  operating  point  coin¬ 
cides  with  the  fuel  cell’s  MPP.  In  order  to  operate  the  fuel  cell 
at  its  MPP  for  every  instant,  the  MPPT  controller  is  made  up  of 
c<P  an  adaptive  feedback  extremum  seeking  algorithm  and  a  switch 
^  mode  power  regulator.  The  purpose  of  the  extremum  seeking 

E  method  is  to  iteratively  adjust  the  fuel  cell’s  current  and  thus 

&  steer  the  operating  point  to  the  maximum  of  the  I-P  plane  which 
corresponds  to  a  maximum  output  power.  The  power  regulator 
forces  the  fuel  cell  to  work  with  the  reference  current  given  by 
the  extremum  seeking  algorithm. 

Extremum  seeking  is  an  adaptive  non-linear  control  method 
which  has  been  used  since  1950s,  but  theoretical  foundations 
for  its  stability  and  performance  were  established  just  recently 
by  Krstic  and  Wang  [10].  The  approach  enjoys  two  main  advan¬ 
tages  when  applied  to  fuel  cell  MPPT  control.  Firstly,  it  does  not 
require  any  prior  knowledge  of  the  fuel  cell;  secondly,  according 
to  the  adaptive  control  law,  the  MPP  will  definitely  be  achieved 
when  the  control  is  convergent. 

The  paper  is  organized  as  follows.  A  dynamic  model  of  PEM 
fuel  cells  is  presented  in  Section  2;  and  in  Section  3  an  adaptive 
MPPT  controller  is  designed.  In  Section  4  extremum  seeking  is 
applied  to  the  fuel  cell  power  system  and  simulation  results  are 
discussed  in  detail. 

2.  Fuel  cell  power  performance  curve  analysis 

2.7.  Fuel  cell  simulation  model 

2.1.1.  Dynamic  gas  transport  model 

The  proportional  relationship  between  the  flow  of  gas  through 
a  valve  and  the  partial  pressure  can  be  stated  as  [1 1] 

<?h2 
Ph2 

and 

<?o2 
Po2 

where  qu2  is  molar  flow  of  hydrogen  (kmol  S-1),  Ph2  hydrogen 
partial  pressure  (atm),  ku2  hydrogen  valve  molar  constant 
(kmol(atmS)-1),  kan  anode  valve  constant  (y/k  mol  kg 
(atms)-1),  Mh2 molar  mass  of  hydrogen  (kg  kmol-1). 

For  hydrogen,  the  derivative  of  the  partial  pressure  can  be 
calculated  by  using  the  following  perfect  gas  equation 

d  RF 

S^2  =  — (3) 

where  R  is  the  universal  gas  constant  ((1  atm)  (kmolK)-1)),  T 
absolute  temperature  (K),  Van  anode  volume  (1),  </h2  hydrogen 
input  flow  (kmolS-1),  q ^  hydrogen  output  flow  (kmolS-1), 
qrU2  hydrogen  flow  that  reacts  (kmol  S-1). 

The  amount  of  hydrogen  consumed  in  the  reaction  can  be 
calculated  from  the  following  electrochemical  principle: 

NI 

<*k  =  ^  =  I  (4) 


a/MH2 


=  ku2 


\j  mo2 


=  ko2 


(1) 


(2) 
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where  N  is  the  number  of  the  series-wound  fuel  cells  in  the 
stack,  I  the  stack  current  (A),  F  Faraday’s  constant  (C  kmol-1), 
kv  modeling  constant  (kmol(sA)-1). 

When  the  output  flow  is  replaced  by  Eq.  (1)  and  the  Laplace 
transform  is  applied  to  Eqs.  (3)  and  (4),  hydrogen  partial  pressure 
can  be  rewritten  in  the  s  domain  as: 


Pu2  = 


l/fa2 

1  +  Tu2S 


(<  -  2 kj) 


where 

_  Van 
™2  kH2RT 


(5) 

(6) 


By  the  same  means,  the  equations  for  oxygen  partial  pressure 
can  be  derived  as: 


Po2  = 


l/*02 

1  +  T0  2S 


(<?S2  ~  kv1) 


where 


(7) 


entropy  change,  the  expression  is  [13]: 

£nemst  =  1-229  -  8.5  x  1( )~\t  -  298.15)  +  4.308 

x  10~57\ln  Ph2  +  0.5  •  In  Po2)  (10) 


Activation  overvoltage  ?7act  is  described  by  the  Tafel  equation, 
which  can  be  expressed  as  [13] 


?7act  —  §l  +  §2  7"  +  §3  7"  In  Cq2  +  §4?"  In  7  (11) 


where  §(/  =  1-4)  are  parametric  coefficients  for  each  cell  model. 
Co2  is  the  concentration  of  dissolved  oxygen  at  the  gas/liquid 
interface  (mol  cm-3),  which  can  be  calculated  by  means  of 


Co2  = 


_ Pfh _ 

(5.08  x  106)  x  exp(— 498/7") 


(12) 


Ohmic  overvoltage  ^0hmic  results  from  the  resistance  of  the 
polymer  membrane  in  electron  and  proton  transfers.  It  can  be 
expressed  as 


ro2  — 


ko2RT 


(8) 


Eq.  (5)  describes  the  relationship  between  stack  current  and 
hydrogen  partial  pressure,  and  Eq.  (7)  the  relationship  between 
stack  current  and  oxygen  partial  pressure.  As  the  load  draws  cur¬ 
rent,  the  reactants  -  hydrogen  and  oxygen  -  become  depleted 
in  the  fuel  cell  stack,  and  both  partial  pressures  drop  accord¬ 
ingly.  To  protect  the  fuel  cell  plant  from  reactants  starvation, 
commonly  excessive  amounts  of  hydrogen  and  oxygen  are  pro¬ 
vided  for  the  stack  (g^  >  2kr7,  <?q2  >  kvI).  A  higher  excess 
ratio  leads  to  higher  partial  pressures,  and  then  a  higher  fuel 
cell  voltage.  However,  too  much  excess  flow  is  a  problem  as 
it  dries  out  the  membrane  and  consumes  much  more  para¬ 
sitic  power.  Pukrushpan  et  al.  [12]  did  a  high  quality  study 
on  how  to  control  the  excess  ratio.  They  designed  several  air 
flow  controllers  to  regulate  the  input  flow  rate  so  that  it  is 
always  twice  as  much  as  the  reaction  rate.  To  focus  on  the 
MPPT  control  problem,  this  paper  will  not  go  into  details 
on  the  excess  ration  control  issue  but  simply  supply  the  fuel 
cell  with  a  constant  flow  rate  that  is  sufficient  for  consump¬ 
tion. 


2.7.2.  Polarization  curve  model 

The  fuel  cell  voltage  as  a  function  of  current  density  in  a 
steady  state  can  be  represented  by  a  polarization  curve,  which  is 
influenced  by  such  parameters  as  the  cell  temperature,  oxygen 
partial  pressure,  hydrogen  partial  pressure  and  membrane  water 
content.  When  current  is  drawn  from  a  fuel  cell,  the  cell  voltage 
Vcell  decreases  from  its  equilibrium  thermodynamic  potential 
7fnernst  (open  circuit  voltage).  This  voltage  drop  consists  of  acti¬ 
vation  loss  ?7act>  ohmic  loss  ?70hmic  and  concentration  loss  ^con- 
The  basic  expression  for  the  cell  voltage  is: 


Ecell  —  T^nernst  +  ?7act  +  ^ohmic  T  Ocon  (9) 

Reversible  thermodynamic  potential  7fnemst  is  described  by 
the  Nernst  equation.  With  literature  values  for  the  standard- state 


V  ohmic  —  I  Pm 


(13) 


The  ohmic  resistance  Rm  is  given  by 

r  r n  try\ 

Pm  =  —7—  (14) 

A 

where  rm  is  membrane  resistivity  (Q  cm)  to  proton  conductivity, 
tm  membrane  thickness  (cm),  A  cell  active  area  (cm2).  Mem¬ 
brane  resistivity  depends  strongly  on  membrane  humidity  and 
temperature,  and  can  be  described  by  the  following  empirical 
expression  [14] 

181.6  [1  +  0.03(7/ A)  +  0. 0062(7/303 )2(// A)2-5] 

y  - - - -  (15^ 

m  [A.m-  0.634  -3(7/ A)]  exp  [4.1 8(7’  -  303/7’)] 

where  Xm  is  the  membrane  water  content.  The  membrane  water 
content  is  a  function  of  the  average  water  activity  am: 


0.043  +  17.81<2m  —  39.85 a^  +  36^,  0  <  am  <  1 
14  +  1.4(<2m  —  1),  1  <  cim  <  3 

(16) 


The  average  water  activity  is  related  to  the  anode  water  vapor 
partial  pressure  Pu?an  and  the  cathode  water  vapor  partial  pres¬ 
sure  Pv,  ca: 

1/  V  1  Pv,SL  n  +  Py,CSL 

am  —  ~(^an+^ca)  —  ~  ~  (17) 

Z  Z  sat 

The  saturation  pressure  of  water  7^  can  be  figured  out  with 
the  following  empirical  expression: 

logio^sat  =  -2.1794  +  0.029537’  -  9.1813  x  10"57’2 

+  1.4454  x  10“77’3  (18) 


The  value  of  Xm  varies  between  0  and  14,  which  is  equivalent 
to  the  relative  humidity  of  0%  and  100%.  Under  supersaturated 
conditions,  however,  the  maximum  possible  value  of  A.m  can 
be  as  high  as  23.  In  addition,  Xm  can  also  be  influenced  by 
the  membrane  preparation  procedure,  the  relative  humidity  of 
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Fig.  2.  PEM  fuel  cell  simulation  model. 


the  feed  gas,  the  stoichiometric  ratio  of  the  feed  gas,  and  the 
age  of  the  membrane  [15].  Hence,  in  this  paper,  Am  is  consid¬ 
ered  as  an  adjustable  parameter  with  a  possible  value  between  0 
and  23. 

Concentration  overvoltage  r]c on  results  from  the  concen¬ 
tration  gradient  of  reactants  as  they  are  consumed  in  the 
reaction.  The  equation  for  concentration  overvoltage  is  shown  by 
[16] 

>?con  =  ^ln(3-2_j  (19) 

nF  \  ilA  J 

where  z'l  is  the  limiting  current.  It  denotes  the  maximum  rate  at 
which  a  reactant  can  be  supplied  to  an  electrode. 

Combining  the  dynamic  gas  transport  model  and  the  static 
polar  curve  model  together  gives  birth  to  a  whole  fuel  cell  sim¬ 
ulation  model,  as  shown  in  Fig.  2.  The  model  parameters  are 
listed  in  Table  1. 


Table  1 

Model  parameters 


Parameter 

Value 

Ref. 

F  (C  kmol-1) 

96484600 

R  (J  kmol-1  K) 

8314.47 

N 

35 

[14] 

A  (cm2) 

232 

[14] 

kr  =  N/4F 

9.07  x  10-8 

Ku2  (kmol  s-1  atm) 

4.22  x  10“5 

[11] 

Ko2  (kmol  s-1  atm) 

2.11  x  10-5 

[11] 

th2(s) 

3.37 

[11] 

T02(s) 

6.47 

[11] 

Imem  (cm) 

0.0178 

[14] 

£t 

-0.944 

[14] 

£2 

0.00354 

[14] 

£3 

7.8  x  10-8 

[14] 

£4 

-1.96  x  10“4 

[14] 

iL  (A  cm-2) 

2.0 

Assumed 

2.2.  Power  performance  curve  analysis 

The  output  power  of  the  fuel  cell  stack  can  be  calculated  by 
the  expression 

^stack  =  ^Fcell  I  (20) 

and  the  power  characteristic  of  the  external  load  R l  is 

r’load  =  I2RL  (21) 

As  an  example,  Fig.  3  shows  the  fuel  cell  stack  power  per¬ 
formance  curves  when  the  membrane  water  content  is  12  and 
16  and  the  load  is  0.065  Q  and  0.03  £2.  The  other  working  con¬ 
ditions  are  the  temperature  of  70  °C,  the  oxygen  pressure  of 
1  atm,  and  the  hydrogen  pressure  of  1  atm.  It  can  be  seen  clearly 
that: 

•  The  operating  point  is  at  the  intersection  of  the  fuel  cell  power 
curve  and  the  load  power  curve.  Thus  the  power  delivered 
takes  a  unique  value  based  on  both  the  load  resistance  and  the 
fuel  cell  operation  parameters. 

•  There  exists  a  unique  maximum  power  point,  at  which  the 
fuel  cell  produces  maximum  power. 

When  the  membrane  water  content  is  12,  fuel  cell  produces  a 
maximal  power  of  3812  W  with  a  load  resistance  of  0.065  ohm. 
If  the  load  resistance  is  changed  to  0.03  ohm,  the  working  point 
moves  from  A  to  B,  and  the  power  drops  from  3812  to  3181 W. 
This  means  only  83%  of  the  power  potential  is  used.  If  the  mem¬ 
brane  water  content  increases  from  12  to  16,  the  working  point 
shifts  from  A  to  C,  and  thus  the  power  rises  from  38 12  to  4735  W, 
97%  of  the  possible  maximal  4864  W. 


Fig.  3.  The  PEM  fuel  cell  power  curves  under  different  operation  conditions. 
Oxygen  pressure  and  hydrogen  pressure  are  both  regulated  to  1  atm. 
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3.  Adaptive  MPPT  controller  design 

3.1.  Schematic  of  the  fuel  cell  MPPT  control  system 

Since  the  power  available  from  the  fuel  cell  is  limited,  it  is 
necessary  to  adjust  the  operation  point  to  settle  on  or  trace  the 
maximum  power  points.  The  majority  of  fuel  cell  control  lit¬ 
erature  is  focused  on  regulating  and  tracking  system  states  to 
known  setpoints  or  trajectories.  Thus  the  specific  MPPs  can  be 
chosen  a  priori  by  solving  a  steady-state  optimization  problem 
subject  to  operational  constraints.  In  this  way,  the  optimality  of 
such  operating  policies  is  completely  dependent  on  the  quality 
of  the  fuel  cell  model.  Besides  the  model  mismatch,  the  fuel  cell 
operation  strongly  depends  in  a  nonlinear  way  on  many  operat¬ 
ing  parameters,  some  of  which  are  even  unmeasurable.  Owing 
to  the  uncertainty  and  time-varying  dynamic  properties  of  the 
electrochemical  processes,  the  MPPs  chosen  a  priori  will  be  sub- 
optimal,  and  can  be  improved  by  performing  the  optimization 
in  real  time.  When  fuel  cell  models  are  not  well  defined  or  when 
only  limited  operation  parameter  measurements  are  available,  it 
is  necessary  to  consider  a  model-free  MPPT  control,  which  is 
able  to  automatically  tune  the  working  point  in  the  right  direc¬ 
tion.  Fig.  4  shows  the  schematic  of  the  proposed  two-level  MPPT 
control  system. 

Changing  the  fuel  cell  operating  parameters  can  shift  the 
working  point,  but  regulating  the  external  load  impedance  is 
more  effective  to  locate  the  MPP.  Thus  the  fuel  cell  operating 
parameters  are  usually  used  to  maintain  good  working  condi¬ 
tions.  An  additional  power  conditioner  between  the  fuel  cell  and 
the  load  can  change  the  equivalent  load  impedance.  A  supervi¬ 
sory  MPPT  controller  functions  primarily  as  an  MPP  searching 
unit,  deciding  the  fuel  cell  working  point.  The  lower-level  DC- 
DC  controller  receives  commands  from  the  supervisory  MPPT 
controller  and  then  generates  detailed  control  instructions  for 
the  power  conditioning  circuit. 


DC-DC  Power 

conditioning  circuit  Power 


Fig.  4.  Schematic  of  the  fuel  cell  MPPT  control  system. 


The  current  mode  power  conditioner  can  efficiently  transform 
a  regulated  DC  output  power  from  a  fluctuating  fuel  cell  source 
power  to  a  variable  load  resistance  in  an  appropriate  form.  The 
converter  with  a  variable  current  gain  G  will  result  in  the  fuel 
cell  seeing  an  equivalent  load  resistance. 

*equ  =  G2RL  (22) 

Therefore  a  feedback  MPP  tracking  algorithm  iteratively  reg¬ 
ulates  the  fuel  cell  current  with  the  power  conditioner,  and  the 
power  conditioner  forces  the  fuel  cell  to  work  at  the  current  by 
shifting  the  equivalent  load  impedance.  When  the  current  coin¬ 
cides  with  the  MPP,  the  equivalent  load  impedance  matches  the 
internal  resistance  of  the  fuel  cell. 

The  current  in  the  DC-DC  converter  is  generally  controlled 
by  a  switching  concept.  The  power  switch  MOSFET  Q  can  be 
turned  on  or  off  in  a  controlled  fashion,  and  the  average  cur¬ 
rent  is  a  function  of  the  on-time  of  the  switch,  the  pulse  width, 
and  the  switching  frequency.  The  most  common  techniques  of 
controlling  a  converter  are  pulse- width  modulation  [17]  and 
sliding-mode  control  [18].  The  power  converter  transient  can 
settle  in  milliseconds  [19],  which  is  sufficiently  brief  to  be  safely 
ignored.  Thus,  the  converter  can  be  treated  as  a  direct  through 
unit  which  is  able  to  instantaneously  regulate  the  fuel  cell  current 
to  commands  from  the  MPP  tracing  algorithm. 

3.2.  Perturbation  and  observe  (P&O)  control  algorithm 

The  fuel  cell  working  point  is  mainly  determined  by  the  super¬ 
visory  MPPT  control  algorithm.  Although  few  discussions  on 
MPPT  control  for  fuel  cells  can  be  found,  a  number  of  MPPT 
methods  have  been  proposed  in  photovoltaic  power  applications. 
Among  them,  the  perturbation  and  observe  (P&O)  algorithm, 
also  known  as  the  “hill  climbing”  method,  is  by  far  the  most 
commonly  used  in  practice  because  of  its  simplicity  in  algorithm 
and  ease  of  implementation.  The  current  based  P&O  algorithm 
tracks  the  peak-power  current  on  the  basis  of  the  past  fuel  cell 
currents  and  powers.  Fig.  5  shows  the  basic  form  of  the  P&O 
algorithm.  The  operating  current  of  the  fuel  cell  is  perturbed  by 
a  small  constant  increment,  and  the  resulting  change  of  power, 
A P,  is  observed.  If  A P  is  positive,  then  it  is  supposed  that  it 
has  moved  the  operating  point  closer  to  the  MPP.  In  this  case, 
further  current  perturbations  in  the  same  direction  will  move 
the  operating  point  toward  the  MPP.  If  A P  is  negative,  then  the 
operating  point  has  moved  away  from  the  MPP,  and  the  direc¬ 
tion  of  perturbation  should  be  reversed  to  move  back  toward  the 
MPP. 

Many  improvements  on  the  P&O  algorithm  have  been  pro¬ 
posed.  For  example,  a  variable  step  size  of  perturbation  with  the 
slope  of  power  as  a  variable  is  as  follows  [20] : 

A  P[k] 

/[k+l]  =  /M  +  M-—^  (23) 

A  I[k] 

where  M  is  the  step-size  corrector,  and  AP[k]/AI[k]  denotes 
the  instantaneous  power  gradient  at  the  Kth  sampling  period.  In 
order  to  find  the  direction  of  change  for  maximizing  the  power, 
the  P&O  method  periodically  perturbs  the  operation  point  of 
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the  system.  If  the  sign  of  power  derivative  A P  and  that  of  volt¬ 
age  derivative  A I  are  the  same,  the  reference  current  should  be 
further  increased,  and  vice  versa. 

A  major  problem  of  the  traditional  numerically  based  MPP 
control  methods  is  that  they  are  lacking  in  solid  theoretical  sup¬ 
port.  As  a  result,  their  stability  and  robustness  is  hard  to  analyze 
and  may  deteriorate  in  some  unfavorable  environments.  The 
P&O  method,  for  example,  requires  the  plant  dynamics  settle 
down  before  optimization.  It  works  well  under  static  conditions 
but  exhibits  erratic  behaviors  under  rapidly  changing  conditions. 
For  instance,  an  abrupt  change  in  membrane  water  content  will 
be  followed  by  an  immediate  change  of  power.  This  leads  to 
AP[k]/AI[k]  ->  oo  and  makes  the  step  size  too  big,  as  shown  in 
Eq.  (23). 

3.3.  Theory  of  the  proposed  adaptive  MPPT  controller  for 
fuel  cell  systems 

Recently,  Krstic  et  al.  [10,21]  presented  a  systematic 
extremum- seeking  control  methodology  supported  by  such 
rigorous  theories  as  averaging  and  singular  perturbation. 
This  real-time  optimization  methodology  involves  a  nonlinear 
dynamic  system  with  feedback  and  adaptation.  These  features 
render  the  extremum-seeking-based  algorithm  more  effective 
than  classical  optimization  and  adaptive  control  algorithms  in 
many  advanced  applications,  such  as  traction  maximization  in 
antilock  braking  for  a  car  [22],  productivity  maximization  of  a 
continuous  stirred  tank  bioreactor  [23],  power  reduction  max¬ 
imization  of  a  flight  [24],  pressure  rise  maximization  of  an 
aeroengine  compressor  [25],  autonomous  vehicle  target  tracking 
[26],  and  PID  tuning  [27]. 

As  compared  with  classical  numerically  based  MPPT 
methods,  the  extremum- seeking  control  approach  enjoys  two 
advantages.  Firstly,  the  optimization  problem  of  power  maxi¬ 
mization  is  explicitly  solved  by  using  the  dynamic  adaptation- 
based  feedback  control  law  with  sinusoidal  perturbation.  The 


MPP  is  hence  guaranteed  to  be  achieved  when  the  control  is 
convergent.  Secondly,  extremum  seeking  is  a  nonmodel-based 
adaptive  control  method  that  can  be  used  to  find  an  unknown, 
optimal  operating  condition  for  the  nonlinear  system,  where  the 
nonlinearity  has  a  local  minimum  or  maximum.  In  other  words, 
our  approach  does  not  require  any  parameterization  or  structural 
formalization  of  the  modeling  uncertainty. 

With  the  self-optimizing  extremum  algorithm  as  the  MPPT 
controller,  the  control  objective  is  for  the  fuel  cell  operating 
point  to  rapidly  trace  the  MPPs  subject  to  uncertainties  and 
disturbances  from  the  fuel  cell  and  the  external  load.  A  block 
diagram  for  extremum  seeking  implemented  on  a  fuel  cell  sys¬ 
tem  is  shown  in  Fig.  6.  7ref  is  our  estimation  of  the  unknown 
MPP  current  7m ax,  at  which  the  fuel  cell  produces  maximum 
power  Pmax.  a ,  co  and  0  are  the  amplitude,  frequency  and  phase 
shift  of  the  probing  signal,  respectively,  co h  is  the  cut-off  fre¬ 
quency  of  the  high-pass  filter,  k  is  the  positive  adaptive  gain  of 
the  integrator.  C(s)  is  the  compensator. 

For  the  sake  of  clarity,  this  paper  only  deals  with  a  controller 
with  C(s)  =  1  and  0  =  0,  as  given  in  Fig.  7.  This  controller  design 
is  the  simplest  form  of  an  extremum- seeking  algorithm  where 
the  algorithm  does  not  employ  any  compensation  or  demodu¬ 
lator  phase  shift.  Although  the  MPPT  controller  appears  to  be 
extremely  simple,  the  analysis  of  the  scheme  is  indeed  intricate. 
The  MPPT  controller  locates  the  MPP  by  injecting  a  small-signal 
periodical  sinusoidal  perturbation  into  the  current  7ref .  If  the  per¬ 
turbation  frequency  is  sufficiently  large,  the  fuel  cell  will  appear 
like  a  static  map  and  its  dynamics  do  not  interfere  with  the  MPPT 
controller.  The  small  perturbation  will  create  a  response  of  P. 
Then  the  high-pass  filter  and  the  integrator  are  applied  to  elim¬ 
inate  the  static  component  of  power  and  extract  the  gradient 
information. 

An  elementary  intuitive  verification  is  given  as  follows.  It 
consists  of  two  steps:  first,  it  is  proved  that  the  MPPT  con¬ 
troller  can  drive  the  fuel  cell  from  any  working  point  to  the 
MPP;  second,  that  the  MPP  is  an  equilibrium  point  is  proved. 


Fig.  6.  MPPT  controller  scheme  for  the  fuel  cell  system. 
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Fig.  7.  Variation  of  the  dP/dl  on  the  P-I  curve. 


P '  >  0.  The  current  increment  kaP'tl 2  is  accordingly  a  positive 
value,  and  the  working  point  will  move  to  MPP.  Contrarily,  if  Iq 
is  bigger  than  /max,  the  current  will  decrease  to  /max  because  of 
kaP'tl  2  <  2. 

To  prove  that  once  the  working  point  arrives  at  the  MPP,  the 
MPPT  controller  can  stabilize  the  system  in  the  small  neighbor¬ 
hood  of  the  MPP,  we  start  with  defining  the  estimation  error  of 
current  as: 

I  e  =  4nax  —  I  ref  (30) 

A  well  designed  MPPT  controller  should  force  e  to  converge 
to  zero.  The  fuel  cell  P(J)  function  can  be  approximated  in  the 
neighborhood  of  the  MPP  by  Taylor  series 

P=Pmm  +  P\I  -  W)  +  |j-(/  -  W)2  +  •  •  •  (31) 

We  also  can  get 

I  /max  =  I  /ref  /<y  —  Cl  sin  COt  Ie  (32) 


The  following  calculations  drive  home  the  point  that  starting 
from  any  point  (Pq  ,  /o)  at  time  to ,  the  MPPT  controller  can  use  the 
small  modulated  perturbation  to  determine  the  slope  of  the  power 
versus  current  curve  and  drive  the  system  to  the  MPP  The  power 
response  to  the  perturbation  signal  at  Po  can  be  approximated 
by  the  Taylor  series 

P" 

P  =  Po  +  P'(/0  +  a  sin  cot)  +  —  (/o  +  a  sin  cot)1  +  •  •  •  (24) 

After  dropping  the  second  and  higher  order  terms,  Eq.  (24) 
can  be  rewritten  as 


P  —  P0  +  P  /o  T  aP  sin  cot  (25) 

The  MPPT  controller  applies  the  high-pass  filter  s/(s+h)  to  the 
output  power  P  to  eliminate  the  slowly  changing  DC  component 
Po  +  ^/o: 

- [P]  —  aP'  sin  cot  (26) 

s  +  h 

This  signal  is  then  demodulated  by  multiplication  with  sin  cot : 

ap>  apf 

£  =  a  P'  sin2  cot  — - cos  2 cot  (27) 

2  2 


With  the  initial  value  of  the  integrator  being  To,  the  estimated 
current  is 


Aef 


§d  t 


kaP' 


/ 


cos  2cotdt 


(28) 


With  Eq.  (32),  Eq  (31)  can  be  written  as 

/  P"  2 

P  =  pmax  +  P  (a  sin  cot  —  Ie)  H - ( a  sin  cot  —  I e)  +  •  •  • 

(33) 

Obviously  P'  =  0  at  the  MPP,  and  by  dropping  the  higher  order 
terms,  Eq.  (31)  can  be  rewritten  as 

/,  =  />max  +  l/,,'(/-/max)2  (34) 

With  Eq.  (32),  Eq.  (34)  can  be  written  as 

P" 

P  =  Pmax  H  ~(a  sin  cot  —  Ie) 

a2P"  P"  2 

—  /max  H - - - 1" 

a2  P" 

—  aP" Ie  sin  cot - cos  loot  (35) 


Estimation  error  Ie  is  a  very  small  value  in  the  neighbor¬ 
hood  of  the  MPP,  thus  its  quadratic  term  I1  in  Eq.  (35)  can 
be  neglected.  The  MPPT  controller  applies  the  high-pass  filter 
s/(s  +  h)  to  the  output  power  P  to  eliminate  the  slowly  changing 
DC  component  Pmax  +  a2P"/ 4: 


- [P]  =  —aP"e  sin  cot  — 

s  +  h 


a2  Pf/ 

- cos  loot 

4 


(36) 


This  signal  is  then  demodulated  by  multiplication  with  sin  cot 


Noting  that  the  second  term  in  Eq.  (28)  is  a  high  frequency 
periodical  signal  around  zero,  which  can  be  attenuated  by  the 
integrator,  we  get 

/ref  =  /(.+  y/yr  (29) 

As  shown  in  Fig.  7,  if  the  value  of  To  is  smaller  than  that  of 
4nax,  the  working  point  is  at  the  left  side  of  the  MPP,  and  thus 


2  p// 


-aP"  1 e  sinz  cot 


azP 


sin  cot  cos  loot 


a2P"  a2  P" 

- Ie  H - Ie  cos  2 cot 


a2  P" 
8 


(sin  cot  —  sin  3  cot) 


(37) 
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The  estimated  current  is 

,  f  ,  kci2P" 
/ref  —  k  I  t  — 


/ 


+ 


/ 


a2P" 


2 

eCOS  2 cot 


edt 


a2  P" 


(sin  cot  —  sin  3 cot) 


d  t 

(38) 


The  second  term  in  Eq.  (38)  can  be  attenuated  by  the  integra¬ 
tor,  thus 

kn2  P"  C 

/ref  = - —  J  Iedt  (39) 

Eq.  (39)  can  be  rewritten  as 
ka-P" 

/ref  = - —4  (40) 


Because  /max  is  a  constant,  we  can  get  /ref  =  —  Ie  from  Eq. 
(30).  Thus 


ka 2  „ 

Ie  =  - P"le 

e  2  * 


(41) 


Since  P(J)  is  a  convex  function  with  maximum,  its  second 
order  derivative  P"  should  be  negative,  and  ka1  P///2<0.  Thus 
we  conclude  that  Ie  0,  which  means  that  the  fuel  cell  current 
converges  to  /max  and  the  MPPT  controller  drives  the  output 
power  to  the  neighborhood  of  Pmax- 

According  to  the  theory  put  forward  by  Krstic  et  al.  [  10] ,  there 
exists  a  ball  of  initial  conditions  around  the  point  (I,P£)  =  (/max, 
Pmax,  0)  and  constants  66 ,  8  and  a  such  that  for  all  co  e  (0,  co), 
8  e  (0,  8)  and  a  e  (0,  a),  the  solution  (I(t),  Pit),  $(t))  exponen¬ 
tially  converges  to  an  O(co  +  8  +  /^-neighborhood  of  that  point. 


3.4.  Design  consideration  of  the  adaptive  fuel  cell  MPPT 
controller 


The  choice  of  the  control  system  parameters  affects  the  sys¬ 
tem  performance.  As  we  can  see  from  the  above  analysis,  the 


overall  fuel  cell  MPPT  control  system  should  have  three  time 
scales: 

•  fastest — the  periodic  perturbation, 

•  medium — the  high-pass  filter, 

•  slow — the  fuel  cell. 

Thus  the  periodical  sinusoidal  perturbation  signal  frequency 
co  and  the  high-pass  filter  cut-off  frequency  co h  are  chosen  as 
follows: 

co  coh  dynamic  speed  of  fuel  cell  with  controller  (42) 

A  higher  co  can  separate  the  perturbation  signal  from  the  fuel 
cell  dynamic  response  clearly.  But  if  co  is  selected  too  high,  the 
convergence  speed  will  become  too  slow.  Thus,  co  should  be 
selected  to  the  degree  that  it  is  just  sufficient  to  satisfy  Eq.  (42). 

The  value  of  Ka 2  has  a  significant  effect  on  the  system  con¬ 
vergence  speed  as  shown  in  Eq.  (41).  A  small  Ka2  slows  down 
the  convergence.  Increasing  the  adaptation  gain  k  or  the  pertur¬ 
bation  signal  amplitude  a  can  improve  the  tracking  speed,  but  the 
system  will  become  more  sensitive  to  noise  or  ripple  contents. 
Therefore  a  moderate  ka 2  should  be  used. 

4.  Simulation  results  and  discussion 

The  complete  fuel  cell  MPPT  control  system  block  diagram 
is  shown  in  Fig.  8.  The  fuel  cell  is  modeled  as  shown  in  Fig.  2, 
while  the  whole  MPPT  control  system  is  written  and  tuned  by 
simulations  on  the  Simulink/Matlab  environment.  The  current 
and  power  of  fuel  cell  stacks  varies  with  their  active  area  and 
number  of  cells,  and  this  variation  will  influence  the  selection  of 
controller  parameters.  To  enlarge  the  scope  of  our  research  appli¬ 
cation,  noting  that  the  current  densities  and  the  power  densities 
of  different  fuel  cells  are  similar,  we  design  the  MPPT  con¬ 
troller  on  the  basis  of  current  density  and  power  density  instead 
of  current  and  power. 

It  is  crucial  to  note  that  all  of  the  static  or  dynamic  fuel  cell 
characteristics  can  be  unknown.  The  fuel  cell  model  presented 


Fig.  8.  Fuel  cell  MPPT  control  system  block  diagram. 
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in  Section  2.1  must  be  viewed  as  a  theoretical  benchmark  to 
illustrate  and  analyze  the  efficiency  of  the  proposed  MPPT  con¬ 
trol  approach.  Even  the  black-box  neuron  network  model  [28] 
or  the  support  vector  machine  model  [29]  can  also  be  used  here. 
It  is  also  to  note  that  a  stably  operated  fuel  cell  plant  needs  a 
group  of  controllers.  In  addition  to  the  excess  ratio  controller 
mentioned  earlier,  water  and  heat  subsystems  also  need  to  be 
properly  controlled.  In  the  following  simulations,  it  is  assumed 
that  they  have  already  been  well  controlled. 

The  basic  simulation  conditions  of  the  fuel  cell  are:  hydro¬ 
gen  input  flow  of  8  x  10-5kmols-1,  oxygen  input  flow  of 
4x  10-5kmols-1,  cell  temperature  of  70  °C  and  membrane 
water  content  of  14.  The  parameters  of  the  MPPT  control  system 
are:  a  =  5,  co  =  2jt,  k  =  14  and  co^  =  tv. 

Fig.  9  shows  the  fuel  cell  MPPT  control  system’s  responses 
to  step  changes  of  the  membrane  water  content.  It  is  important 
to  understand  that  these  abrupt  changes  of  the  membrane  water 
content  are  for  testing  the  dynamic  response  of  the  MPPT  control 
system,  and  do  not  necessarily  represent  changes  in  a  real  world 
case.  The  membrane  water  content  abruptly  changes  from  8  to 
12  at  the  100th  second,  from  12  to  16  at  the  150th  second  and 
from  16  back  to  12  at  the  200th  second,  as  shown  in  Fig.  9a.  The 
corresponding  MPPs  also  move  from  point  A  to  B,  from  B  to  C 
and  then  from  C  back  to  B,  as  shown  in  Fig.  10a. 

Figs.  9  and  10  show  that  the  MPPT  controller  can  steadily 
drive  the  system  to  the  new  MPP  when  an  abrupt  change  of  the 
MPP  occurs.  For  example,  when  the  membrane  water  content 
abruptly  changes  from  8  to  12  at  the  100th  second,  the  fuel  cell 
operating  point  immediately  shifts  to  point  A\  and  the  MPP 
shifts  to  point  B'.  If  the  P&O  algorithm  governed  by  Eq.  (23)  is 
used,  the  step  change  of  power  from  point  A  to  A'  will  lead  to 
AP[k]/AI[k ]  ->  oo,  and  then  the  P&O  algorithm  will  collapse. 
The  MPPT  controller,  however,  can  converge  to  the  maximum 
power  point  B  smoothly,  as  shown  in  Fig.  9f  and  10b. 

The  movement  from  point  B'  to  B  is  caused  by  the  couple 
between  current  and  oxygen/hydrogen  pressures,  as  described 
by  Eqs.  (5)  and  (7).  When  the  current  rises,  the  depletion  of 
the  reactant  gases  increases,  and  the  gas  pressures  decrease,  as 
shown  in  Fig.  9d  and  e.  And  when  the  reactant  pressures  change, 
the  location  of  the  MPP  changes  accordingly,  as  described  in 
Section  2.1.2.  Strictly  speaking,  the  membrane  water  content  and 
the  fuel  cell  temperature  both  vary  with  current.  When  current  is 
changed  to  seek  the  MPP,  the  location  of  the  MPP  changes  cor¬ 
respondingly.  But  the  MPPT  control  system  is  robust  enough  to 
trace  the  locus  of  MPPs  until  the  equilibrium  point  B  is  reached. 

Fig.  11  shows  the  fuel  cell  MPPT  control  system’s  ramp 
response  to  the  temperature.  The  fuel  cell  temperature  rises 
from  20  to  90  °C  at  the  rate  of  0.1  °Cs_1.  When  the  temper¬ 
ature  changes  from  20  to  90  °C,  the  MPP  moves  from  Point  A  to 
B.  The  MPPT  control  system  traces  the  theoretical  MPP  locus 
at  the  same  time.  The  tracing  error  is  mainly  caused  by  the  slow 
change  of  the  gas  pressures. 

Although  this  MPPT  controller  leads  to  a  forced  oscillation 
current,  it  achieves  a  very  low  power  oscillation  around  the  MPP. 
A  periodical  perturbation  current  of  asm  cot  is  injected  to  the 
fuel  cell  current.  This  leads  to  a  ripple  power,  whose  amplitude 
depends  on  the  relative  location  of  the  operation  point  to  the 
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Fig.  9.  Responses  of  the  MPPT  controller  to  the  step  changes  of  membrane  water 
content:  (a)  step  changes  of  membrane  water  content;  (b)  responses  of  current 
density;  (c)  responses  of  cell  voltage;  (d)  responses  of  hydrogen  pressure;  (e) 
responses  of  oxygen  pressure;  (f)  responses  of  power  density. 

MPP.  The  power  oscillation  amplitude  can  be  approximated  by 
d  P 

Power  amplitude  =  x  a  (43) 

d  I 

As  mentioned  earlier,  dP/dl  — >►  0  at  the  MPP.  Thus  the  power 
amplitude  will  converge  to  zero  at  the  MPP.  Fig.  10b  shows  that 
the  power  ripple  becomes  lower  and  lower  when  the  operation 
point  gets  close  to  MPP  B.  Fig.  12  gives  a  clearer  picture  when 
the  operation  point  is  very  close  to  MPP  B.  The  current  ripple 
factor  in  Fig.  12a  is  about  1%,  but  the  corresponding  power 
ripple  factor  is  only  0.02%  as  shown  in  Fig.  12b. 
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Fig.  10.  Power  responses  in  phase  plan  of  power  density  versus  current  density:  (a)  power  density  responses  in  I-P  plan;  (b)  oscillations  nearby  point  B.  Dot  lines 
are  the  polar  curves  at  different  membrane  water  contents. 


Fig.  11.  Response  of  the  MPPT  controller  to  the  ramp  change  of  temperature:  (a)  power  density  responses  in  time  domain;  (b)  power  density  responses  in  I-P  plan. 


Fig.  12.  Ripples  of  Fuel  cell  MPPT  control  system. 


5.  Conclusion 

A  modeling,  control,  and  simulation  study  of  the  fuel  cell 
MPPT  control  system  is  reported  in  this  paper.  Firstly,  a  the¬ 
oretical  analysis  of  the  power  curve  characteristic  of  the  fuel 
cell  is  made.  The  fuel  cell  has  a  unique  MPP,  at  which  it  pro¬ 


duces  maximum  power.  The  location  of  the  MPP  is  decided  by 
fuel  cell  operation  conditions  in  a  nonlinear  way.  Due  to  the 
complexity  of  the  electrochemical  process,  however,  it  is  hard 
to  predict  whether  the  fuel  cell  is  working  at  the  MPP  by  a 
modeling  approach. 

Therefore,  an  adaptive  MPPT  controller  using  the  extremum¬ 
seeking  algorithm  is  designed  to  automatically  keep  the  fuel  cell 
working  at  the  MPP  all  the  time.  The  controller  is  based  on  solid 
dynamic  adaptive  control  theories.  The  MPP  is  guaranteed  to 
be  achieved  when  the  control  is  convergent.  Also,  this  MPPT 
controller  is  a  model-free  method  which  does  not  require  any 
prior  knowledge  of  the  fuel  cell.  The  validity  of  the  proposed 
fuel  cell  MPPT  control  scheme  over  a  large  operating  range  was 
verified  through  simulations  in  the  MATLAB/Simulink  envi¬ 
ronment.  Simulation  results  show  that  the  system  can  respond 
rapidly  to  abrupt  changes  of  operation  conditions. 

The  proposed  control  scheme  can  not  only  trace  the  MPP  but 
also  online  solve  other  fuel  cell  optimization  problems,  such  as 
maximize  efficiency  or  combination  of  various  objectives.  We 
expect  that  it  can  be  used  in  a  large  scope. 
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